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First measurements of the decays of the three Xcj states to ppK^K~ final states are presented. 
Intermediate — > K^K^ and A(1520) pK^ resonance states are observed, and branching 
fractions for Xcj — 5- pA"''A(1520), A(1520)A(1520), and (j)pp are reported. We also measure branching 
fractions for direct Xcj — >■ ppK^ K~ decays. These are first observations of Xcj decays to unstable 
baryon resonances and provide useful information about the XcJ states. The experiment uses samples 
of XcJ mesons produced via radiative transitions from 106 million ^' mesons collected in the BESIII 
detector at the BEPCII e^e~ collider. 

PACS numbers: 13.25.Gv. 14.20.Pt. 14.40.Be 
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I. INTRODUCTION 

Experimental studies on charmonia decay properties are useful for testing perturbative QCD models and QCD- 
based calculations. In the standard quark model, the Xcj{ J = 0,1,2) mesons are P-wave quarkonium states with 
spin parity 0'*"'', 1'*"'' and 2"*'+. Although they cannot be produced directly in e+e" collisions, radiative decays of the 
"0' into each Xcj occur about 9% [1] of the time and provide large Xcj samples that have proven to be a very clean 
environment for studies of the XcJ states. 

The color octet mechanism (COM) has been shown to play an important role in describing these P-wave quarkonium 
decays [2-5]. Many COM predictions for Xcj decays into meson pairs and pp pairs are in agreement with earlier 
experimental results. However, the predictions for some baryon-antibaryon decays disagree with measured values, in 
particular Xcj ^ AA [6]. At present, only ground state baryons have been observed in Xcj decays [1]. To further 
test COM predictions for P-wave charmonia decay, measurements of excited baryon pair decays are important. This 
paper presents a study of Xcj hadronic decays and measurements of Xcj A(1520)A(1520) decaying to ppK^K~ , 
based on 106 million i}}' events collected with BESIII at BEPCII. The observation of such excited baryon production 
can provide constraints on models of P-wave charmonia hadronic decay. 

II. DETECTOR 

BEPCII [7] is a double-ring e+e~ collider designed to provide a peak luminosity of 10^'^ cm~^s~^ at the center 
of mass energy of 3770 MeV. The BESIII [7] detector has a geometrical acceptance of 93% of Air and has four main 
components: (1) A small-cell, helium-based (40% He, 60% CsHg) main drift chamber (MDC) with 43 layers providing 
an average single-hit resolution of 135 /im, and charged-particle momentum resolution in a 1 T magnetic field of 0.5% 
at 1 GeV/c. (2) An electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) crystals in a cylindrical structure 
(barrel) and two endcaps. The energy resolution at 1.0 GeV/c is 2.5% (5%) in the barrel (endcaps), and the position 
resolution is 6 mm (9 mm) in the barrel (endcaps). (3) Particle Identification is provided by a time-of-flight system 
(TOF) constructed of 5-cm-thick plastic scintillators, with 176 detectors of 2.4 m length in two layers in the barrel 
and 96 fan-shaped detectors in the endcaps. The barrel (endcap) time resolution of 80 ps (110 ps) provides 2cr K/tt 
separation for momenta up to ^ 1.0 GeV/c. (4) The muon system (MUC) consists of 1000 m^ of Resistive Plate 
Chambers (RPCs) in nine barrel and eight endcap layers and provides 2 cm position resolution. 
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III. MONTE-CARLO SIMULATION 

Montc-Carlo (MC) simulation of the full detector is used to determine the detection effieiency of each channel, 
optimize event selection criteria, and estimate backgrounds. The simulation program, BOOST, provides an event 
generator, contains the detector geometry description, and simulates the detector response and signal digitization. 
Charmonium resonances, such as the -0', are generated by KKMC [8, 9], which accounts for effects such as initial state 
radiation and beam energy spread. The subsequent charmonium meson decays are produced with BesEvtGen [10, 11]. 
The detector geometry and material description and the tracking of the decay particles through the detector including 
interactions are handled by Geant4. 

IV. EVENT SELECTION 

Charged tracks must have their point of closest approach to the beamline within ±10 cm of the interaction point 
in the beam direction and within 1 cm of the beamline in the plane perpendicular to the beam and must have the 
polar angle satisfy | cos 6*1 < 0.93. The TOF and energy loss dE/dx measurements are combined to calculate particle 
identification (PID) probabilities for pion, kaon, and proton/anti-proton hypotheses, and each track is assigned a 
particle type corresponding to the hypothesis with the highest confidence level (C.L.). Finally, four tracks identified 
as p, p K^, and are required. 

Photon candidates are selected by requiring a minimum energy deposition of 80 MeV in the EMC. EMC cluster 
timing requirements suppress electronic noise and energy deposits unrelated to the event. 

Kinematic fitting that utilizes momentum and energy conservation is applied under the hypothesis "0' — ?> 7Xc,7 ~^ 
jppK^ K~ . For events with more than one photon candidate, the combination with the smallest xic considered 
for further analysis. 

V. DATA ANALYSIS 

After candidate event selection, distinct Xcj signals are observed in the ppK^ invariant mass distribution, as 
shown in Fig. 1(a). By combining final state particles {p,p, , K^), the A(1520); A(1520) and <j> intermediate states 
can also be seen in the pK^, pK^ and K^K^ invariant mass distributions, as shown in Fig. l(b,c,d), respectively. 
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FIG. 1: Invariant mass distributions of (a) ppK^K , (b) pK , (c) pK^ and (d) K^K 



A. Background Studies 



A sample of 100 million inclusive ip' MC events is used to investigate possible backgrounds. No background events 
survive after candidate selection. 

Potential physics background contributions due to undetected or fake photons and particle misidentification can 
come from the processes: ip' ^ ir^ppK^ , ip' IXcj ~> , jK^K^tt^tt^ , ^ppK^ir^ and ip' — > 

ppK^K" . We produced 2 x 10^ MC events for the first process and 1 x 10^ MC events for each of the other processes 
in order to study these backgrounds. After applying the event selection criteria to MC events, 265 events survive, 
and all of them are from ip' — >■ n^ppK^ . Since the branching fraction of this channel has not been reported by 
the PDG [1], we determine it from our data sample and use the result to estimate the background contribution to be 
about 1.4 events. In addition, a 42.9 pb~^ data sample collected at 3.65 GeV is used to investigate possible continuum 
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backgrounds, and no events survive candidate selection. 



B. XcJ ppK+K- 

The branching fractions for Xcj ppK^K~ are measured excluding the evident A(1520), A(1520) and cj) intermedi- 
ate states seen in Fig. l(b,c,d), by vetoing these events with the mass requirements \M{pK^) — 1.52| > 0.07 GcV/c^, 
\M{pK+) - 1.52| > 0.07 GeV/c2 and \M{K+K-) - 1.02| > 0.03 GeV/c^. The fit of the invariant mass distribution 
of the remaining candidate events is shown in Fig. 2. 
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FIG. 2: The ppK+K' invariant mass distribution and fit result after excluding A(1520), A(1520) and (p intermediate states, 
where dots with error bars are data and the solid curves show the fit result. The dashed line, barely perceptible, is the estimated 
background component. 



The ppK^ mass distribution is fitted with Breit-Wigner functions convolved with Gaussian resolution functions 
to describe the XcJ signals and a flat distribution for the background, as shown in Fig. 2. The Breit-Wigner parameters 
and Gaussian instrumental resolutions are floated in the fit, with the XcJ widths fixed according to PDG [1] values. 
The instrumental resolutions are found to be about 4 MeV/c^, and the masses of the Xc.j from the fit are consistent 
with PDG values within Icr. The observed numbers of events, denoted as Nobs, for XcJ ^ ppK^K~ decays are listed 
in Table I. The branching fractions are calculated according to : 



Bixcj^PpK+K-) 



N^, ■ B{^P' ^ ^Xcj) ■ e' 

where 7V^/ is the total number of tp' events, which is measured to be 106 x 10^ with an uncertainty of 4% [12], the 
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V-' ^ IXcJ branching fractions are taken from PDG [1] to be (9.62 ± 0.31)%, (9.2 ± 0.4)% and (8.74 ± 0.35)% for 
XcOj Xcl and Xc2, respectively, and the detection efficiencies, e, are determined individually for simulated ^/j' decays 
to ppK^K^ via the XcOi Xci and Xc2 states. The results are summarized in Table I. 

TABLE I: The branching fractions for Xcj ppK^ , where errors are statistical only. 
Quantity XcO Xa Xc2 

Nabs 48.2 ±7.7 81.5 ±9.2 131 ± 12 

e(%) 3.8 ±0.1 6.2 ±0.1 6.8 ±0.1 

Bixcj ^ PpK+K-) (10-*) 1.24 ±0.20 1.35 ± 0.15 2.08 ±0.19 



C. XcJ ~^ pAT"*" A(1520) + charge-conjugate (c.c.) 

For analysis of the intermediate states described below, the following invariant mass selection criteria are imposed: 

Xco ■■ 3.365 GeV/c2 < M{ppK+K-) < 3.455 GeV/c^ 
Xcl : 3.490 GeV/c2 < M{ppK+K-) < 3.530 GeV/c^ 
Xc2 : 3.530 GeV/c2 < M{ppK+K-) < 3.580 GeV/c^ 

The three-body decay branching fractions Xcj — > ?3i^^A(1520) + c.c. are measured after rejecting the A(1520) 
for pA'+A(1520) with the requirement \M{pK+) - 1.52| > 0.07 GeV/c^ or for the charge- conjugate A(1520) to 
pK~K{\^2Q) by \A4{j)K~) — 1.52| > 0.07 GcV/c^. The fitted pK^ ± c.c. invariant mass distributions are shown in 
Fig. 3. The fits use Breit-Wigner functions convolved with Gaussians for the signals and Chebyshev polynomials for 
backgrounds, where the Breit-Wigner masses and Gaussian instrumental resolution are free parameters, and the widths 
of the resonances are fixed to their PDG [1] values. The observed numbers of events. Nobs, for Xcj p/\ ^A(1520)±c.c. 
arc shown in Table II. 

The branching fractions are calculated according to: 

^ PA-A,1520) + „J . „,.BW.^,,./gwi520)^pA-).. - 

where the detection efficiencies e are determined by MC simulation of decays to p/v +A(1520) ± c.c. for each of the 
XcOj Xcl and Xc2 states. The results are summarized in Table II. 
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FIG. 3: Invariant mass distributions and fits to pK~ + c.c. in the decays of (a) XcO, (b) Xci and (c) Xc2. Dots with error bars 
are data. Solid lines are results of the fit, and dashed curves represent the background. 



TABLE II: The branching fractions for Xcj — > p_ft'^A(1520) + c.c, where errors are statistical only. 



Quantity 



Xc2 



Nobs 

e{%) 

S(A(1520) -^pK-){%) 



62 ± 12 48 ± 10 79 ± 13 
9.0 ±0.1 12.1 ±0.1 12.4 ±0.1 
22.5 22.5 22.5 



^(XcJ ^ P^i'+A(1520) ± c.c.) (10"") 3.00 ±0.58 1.81 ±0.38 3.06 ± 0.50 



D. XcJ A(1520)A(1520) 

A scatter plot of the invariant mass M{pK~^) versus M{pK~) is shown in Fig. 4(a), where a signal for XcJ 
A(1520)A(1520) is evident. Events remaining after rejecting (ppp events with the veto requirement \AI{K^K^) — 
1.02| > 0.03 GeV/c2 and satisfying \M{pK+) - 1.520| < 0.05 GcV/c^ and \M{pK-) - 1.520| < 0.05 GeV/c^ are 
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selected as candidate A(1520)A(1520) events. Two-dimensional mass sideband regions used to estimate background 
in the signal region S arc indicated by the regions A,B,and C in Fig. 4(a). The events in these sideband regions for data 
are scaled by factors that are determined by the ratios of events in the signal region S to those in the sideband regions 
for MC samples for the background channels, namely -0' — > 7Xc,/ 7p/^''"A(1520), V'' ~^ IXd ~^ 7pA"~A(1520) and 
■0' — 7Xcj ~> ^ppK~^ K~ , are to estimate the background in the signal region S of the data. 
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FIG. 4: (a) Scatter plot of M(pK^) versus M{pK~); (b) Invariant mass spectrum and fits to ppK^ K~ , where dots with error 
bars are events from the signal region. The solid line is the fitting curve for the events from signal region, and the dashed lines 
represent background estimated from the two-dimensional mass sidebands of regions "A, B, C" as shown in (a). 

The mass spectra obtained from the signal and scaled sideband background events in Fig. 4(a) are simultaneously 
fit using Brcit-Wigncr functions convolved with Gaussian resolution functions. The Breit-Wigner masses and the 
instrumental resolutions used for the Gaussians are left as free parameters in the fit. Other background is described 
by a flat distribution. The differences between the results of the fits to the signal and scaled sideband events, shown 
in Fig. 4(b), are used to extract the Xcj -> A(1520)A(1520) yield. We find 28.1 ±9.8 events for Xco ^ A(1520)A(1520) 
and 28.9 ± 7.4 events for Xc2 ^ A(1520)A(1520). No distinct Xci A(1520)A(1520) signal is observed, and a 90% 
C.L. upper limit is given using the Bayesian method. 

The branching fractions are calculated according to: 
S(xcJ ^ A(1520)A(1520)) 



N^, ■ B{^' -> 7Xcj) ■ i3(A(1520) ^ pR-) ■ S(A(1520) pK+) ■ e 
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and the upper limit at the 90% C.L. is calculated as: 
B{xci A(1520)A(1520)) 

Nobs 

N.^, ■ B{^' 7Xci) ■ B(A(1520) -> pR-) ■ B(A(1520) ^ pK+) • e ■ (1 - a^ys) ' 
where the detection efficiencies are determined from MC simulation, which assumes an angular distribution of 1 + 
acos^ 9 for the two-body decays, and the value for a is estimated by fitting the cos 6' distribution of data separately 
for the XcOi Xci and Xc2 states, 9 is the polar angle of a particle in the rest frame of its mother particle, and asys 
denotes the systematic error (discussed below). The results are summarized in Table III. 

TABLE III: The branching fractions for XaJ ~^ A(1520)A(1520). The errors are statistical only, and the upper limit is at the 

90% C.L. 

Quantity XcO Xci Xc2 

Nobs 28.1 ±9.8 < 6.9 28.9 ± 7.4 

e(%) 17.1 ±0.1 16.3 ±0.1 12.2 ±0.1 

S(A(1520) pK){%) 22.5 22.5 22.5 

B(Xcj ^A(1520)A(1520)) (10"'') 3.18 ±1.11 < 0.86 5.05 ± 1.29 



E. x<^j pp4> 

The K^K~ invariant mass distributions and fits to the spectra are presented in Fig. 5 for the XcOi Xcii and 
Xc2- (j) signals are observed clearly in the decays of Xco (Fig. 5(a)) and Xc2 (Fig. 5(c)). The fits use Brcit-Wigner 
functions convolved with Gaussians for the signals, where the Breit-Wigner masses and instrumental resolutions 
are free parameters and resonance widths are fixed at their PDG [1] values, and a Chebyshev polynomial for the 
background. The observed numbers of events. Nobs, for Xcj ~> PP4' are listed in Table IV, as well as an upper limit 
at the 90% C.L. for Xci PP^* using a Bayesian method. The branching fractions are estimated as: 

AT . 

BiXcj PP(f)) 



N^, ■ B{^' ^ 7Xcj) • B((^ ^ K+K-) 
and the upper limit at the 90% C.L. is calculated as: 

Nobs 



S(xci -^pp4>) < 



N^, ■ B{i'' ^ 7Xci) ■ B{4> ^ K+K-) •£•(!- a^ys) ' 
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where detection efficiencies are determined from MC simulation as described above. The results are summarized in 
Table IV. 




M(K" K") (GeV/c^) 

FIG. 5: K~ invariant mass distributions and fits to the spectra in the decays of (a) XcO and (c) Xc2. Dots with error bars 
are data, and solid lines represent the fit results. Dashed curves are background shapes. For decays of (b) Xci, <i> is not seen 
clearly, and the upper limit at the 90% C.L. is given. 

TABLE IV: The branching fractions for Xd ^ PP4'- The errors are statistical only, and the upper limit is at the 90% C.L.. 

Quantity XcO Xci Xc2 

Nobs 42.4 ±8.2 < 13.3 24.4 ±6.8 

e(%) 13.9 ±0.1 17.7 ±0.1 17.7 ±0.1 

B{(l> ^ K+K'){%) 48.9 48.9 48.9 

S(xcJ ->PP9i) (10^^) 6.12 ± 1.18 < 1.58 3.04 ± 0.85 
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VI. SYSTEMATIC UNCERTAINTY 

The main contributions to the branching fraction systematic uncertainties originate primarily from the tracking, 
particle identification, photon reconstruction, kinematic fit, branching fractions of the intermediate states (from 
PDG [1]), total number of V'' events, fitting procedure and the event generator. The contributions of each item are 
summarized in Table V for Xcj ^ ppK+R-, A(1520)A(1520) and Table VI for pA'+A(1520) + c.c. and ppcj). 

From analyses of very clean J/tJj — > K*K and J /^j — !■ pp7r+7r~ decays, the tracking efficiency for MC simulated 
events is found to agree with that determined using data to within 2% for each charged track. Hence, 8% is taken as 
the systematic uncertainty for the four charged track final state. 

The candidates of the selected final state require tracks be identified as p, p, or . Comparing data and MC 
event samples for J /tjj — > ir^TT^pp and J /tp — > K* K , a difFerence in MC and data particle identification efficiency of 
2% is obtained for each particle. Hence, 8% is taken as the systematic uncertainty for ppK^K^ identification. 

Photon reconstruction efficiency is studied using -0' tt+tt" J/'0 — > 77r+7r~pp, and the difFerence between data and 
MC is about 1% per photon [12]. 

To estimate the uncertainty from kinematic fitting, a -0' — >■ 7Xcj ^ 7pp7r+7r~ sample is selected to study efficiency 
differences between data and MC. Errors of 1.4%, 1.6% and 2.3% are obtained for decays of XcO> Xci and Xc2, 
respectively. 

Uncertainties due to the decay model used in simulation for two-body and three-body decay channels are estimated 
by varying the a values in the decay angular distributions 1 + a cos^ 6. For two-body decay channels, a is varied over 
a range such that the angular distribution in MC is consistent with that of data. For three-body decays, the accuracy 
of the angular distributions in data are limited by low statistics. To be conservative, we vary a from —1 to 1 and the 
resulting differences are taken as the systematic uncertainty. 

Uncertainties in the fitting procedure are obtained by altering background shapes and fit intervals. Uncertainties 
from the mass window requirements, obtained by changing the requirements, of Xcj, A(1520), A(1520) and (p are 
shown. 

Uncertainties in the reconstruction efficiency for XcJ ppK^K^ due to other possible intermediate states, Xci — >■ 
p_ft'+A(1600) -I- c.c. and XcOi Xc2 — ^ pK^ h^{\67Q) + c.c, which are not pronounced in the data, are summarized in 
Table V. Both masses and widths of A(1600) and A(1670) are poorly determined, and their branching fractions are 
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not available. Their branching fractions are taken conservatively as 5 x 10~^, and the systematic uncertainties are 
the differences between with and without the intermediate states. 

The total number of ip' events with an uncertainty of 4% is obtained by studying inclusive hadronic ip' decays [12]. 
The total systematic uncertainty is obtained by summing up uncertainties contributed from all individual sources in 
quadrature. 

TABLE V: Systematic uncertainties expressed in percent (%) for the decay modes Xcj ppK+K~ and Xd ^ A(1520)A(1520). 





XcJ - 


ppK+R- 




XcJ - 


^ A(1520)A(1520) 






XcO 




Xc-2 


XcO 


Xcl 


Xc2 


Tracking 


8.0 


8.0 


8.0 


8.0 


8.0 


8.0 


PID 


8.0 


8.0 


8.0 


8.0 


8.0 


8.0 


Photon recon. 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Kinematic Fit 


1.4 


1.6 


2.3 


1.4 


1.6 


2.3 


Fitting 


1.5 


0.6 


0.0 


0.0 


0.0 


0.0 


Mass window 


4.2 


6.0 


2.8 


8.2 




11.0 


a value 








3.3 




4.1 


Branching fraction 


3.2 


4.3 


4.0 


7.0 


7.6 


7.4 


N^' 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Efficiency 


5.0 


1.2 


5.9 








Total 


14.2 


14.3 


14.5 


16.6 


14.3 


18.5 



VII. RESULTS AND DISCUSSION 

The measured branching fractions for the twelve decay modes decaying to ppK^ K~ are summarized in Table VII. 
From the 106 million -0' decays observed by BESIII at BEPCII, we report first measurements of these branching 
fractions with uncertainties ranging from 20% to 40%. With larger statistics in future BESIII running, we ex- 
pect to improve these measurements and to be able to observe A(1520)A(1520) in Xci decays. The excited baryon 
A(1520)A(1520) decays provide new information for evaluating model predictions of Xcj hadronic decays. 
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TABLE VI: Systematic uncertainties expressed in percent (%) for the decay modes Xcj — >■ pK^ A{1520) + c.c. and Xcj — >■ PP<^- 





XcJ - 


pA A(io20j + 


c.c. 




XcJ PP(t> 






XcO 


Xcl 


Xc2 




Xcl 


Xc2 


Tracking 


8.0 


8.0 


8.0 


8.0 


8.0 


8.0 


PID 


8.0 


8.0 


8.0 


8.0 


8.0 


8.0 


Photon recon. 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Kinematic Fit 


1.4 


1.6 


2.3 


1.4 


1.6 


2.3 


Fitting 


9.4 


5.9 


6.8 


4.5 




4.7 


Mass window 


2.2 


3.6 


8.8 


2.1 




1.0 


In o 






2.2 


4 n 


^ Q 


9 


Branching fraction 


5.4 


6.2 


5.9 


3.4 


4.4 


4.1 




4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Total 


16.6 


15.5 


17.7 


14.1 


13.5 


14.0 



TABLE VII: Summary of branching fractions for twelve Xcj decay modes to ppK^K . The first errors are statistical, and the 
second ones are systematic. The upper limits are at the 90% C.L. including the systematic errors. 

XcO Xcl Xc2 

B(xcj ^ PpK+K-) {10-'^) 1.24 ±0.20 ±0.18 1.35 ± 0.15 ± 0.19 2.08 ± 0.19 ± 0.30 

S(XcJ ^P-R'+A(1520) ± c.c.) (10-") 3.00 ±0.58 ±0.50 1.81 ± 0.38 ± 0.28 3.06 ± 0.50 ± 0.54 

S(Xc.7 A(1520)A(1520)) (10-*) 3.18 ± 1.11 ± 0.53 < 1.00 5.05 ± 1.29 ± 0.93 

BiXcJ ^ PP(I>) (W-^) 6.12 ±1.18 ±0.86 < 1.82 3.04 ± 0.85 ± 0.43 
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